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Abstract 





The Lyman and Balmer emission of a planetary corona depend on the 
exospheric temperature, the inteqrated column density of solar-illuminated 
hydrogen, and the region of phase space occupied by particles. Measure- 
ments of the intensity alone are incapable of defining the exosphere 
unambiguously. Line profiles, with high spectral resolution, can show 
whether a non- thermal component of the escaping hydrogen is present and 
can indicate at what altitude satellite orbits of hydrogen atoms are 
depleted. It is necessary, however, to plan the observations carefully 
if they are to be fitted usefully to a model. 


I. Introduction 

Two matters of substance complicate the interpretation of direct 
optical observations of a planetary coronas the component temperatures 
t the degree to which satellite orbits are present are generally 
unknown . 

The temperature problem evidently exists for the earth's corona, 
even though independent scale-height measures from satellite drag are 
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available. The difficulty is that intensity observations of Lya require a 
higher escaping flux of hydrogen than can be accounted for by thermal escape 
of H atoms (Hunten, 1973 ? Hunten and Strobel, 1974 ? Liu and Donahue, 1974; 

Bertaux, 1975; Tinsley, et al., 1975; Thomas and Anderson, 1975). [More pre- 
cisely, observed variations of the temperature and intensity of Lya imply 
large time variations in the thermal evaporation of H atoms. But the total 
loss^hydroger must be constant, because the upward flow is limited by eddy 
diffusion much lower in the atmosphere. Hence additional escape that compensates 
for the fluctuations in thermal evaporation is implied.] It becomes necessary 
to invoke an additional mechanism, either charge exchange of cold H with hot H 
{Cole, 1966; Tinsley, 1973; Tinsley et al., 1975) or a chemical mechanism that 
releases H atoms with high velocity. To the extent that the non-thermal mechanism 
is important, the ordinary exospheric scale height is unhelpful. 

The temperature problem may also exist for Venus (McElroy, 1969? Kumar 
and Hunten, 1974) , where Mariner IV observations suggested a superposition of 
two intermingled coronae with distinct temperatures. 

The problem of satellite orbits is equally unsettling, as we have neither 
observation nor theory to indicate the extent that such orbits are populated. 

My earlier paper (Chamberlain, 1963, Sec. 5b) suggested an equilibrium might 
be established by the creation of satellites through collisions and their destructic 
through both collisions and photoionization. More likely, radiation pressure, 
which can distort circular orbits into elliptical ones and lower their perigees 
(Slowey, 1969; Thomas and Bohlin, 1972? Bertaux and Blamont, 1973), is the dominant 
destruction mechanism. If so, variations of solar Lya could alter the effective 
base level of satellite orbits, so that no single distribution of orbits would 
be universally appropriate. 

Therefore, it seems likely that the few data we have need to be 
supplemented. The measurement of line profiles with high spectral 
resolution might clarify matters. The basic theory for profiles was 
developed earlier (Chamberlain, 1963, Sec. 7) but no numerical calcula- 
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tions were made. Here I wish to indicate the kind of profile that can 
be expected for the simplest possible situation: a single temperature, 
radially outward observations, and a well-defined "satellite critical 
level" (the maximum perigee for satellite orbits) . The theory can be 
generalised for a two-component gas, but the examples given here should 
suffice to outline the problem. The only measurements available to date 
(Atreya, Hays, and Nagy, 1975) were not interpreted with allowance for 
the vacant orbits, and the long time-integration of their observations 
would in any case make it difficult to construct an appropriate coronal 
model. In principle, the spectrum could be computed for an arbitrary 
direction of observation, but the problem is enormously simplified if 
the observations are made toward the zenith. 


II. Profile Theory: A Summary 
The radial Doppler profile is 


4ffl(v r , X x ) = g(Lya) 


-X -V 

N c r X e 
c c c 
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1/2 


Stt, \). 


( 1 ) 


Here I is the photon intensity for a Doppler displacement v r » when the 
shadow height intersects the line of sight at distance r^ from the 
planet's center. The dimensionless units used in the theory are 


_ e JE. ZL 

~ kT r ' 5 = U ' V 2 
c U 


( 2 ) 


where U = 2kT c /M, M is the mass of the planet, M the mass of the hydro- 
gen atom, and M and T the H density and the Maxwellian temperature, 
c c 

respectively, at the critical level, r . The g(Lya) -factor is the 

c 
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number of solar photons scattered in all directions per second per H 
atom (Chamberlain 1961, p 425) . Me neglect multiple scattering, assuming 
an optically thin corona. An added refinement would include the anisotropy 
of resonance scattering in the g-factor (Brandt and Chamberlain 1959) . 

The factor SfC,^) in equation (1) is the "radial spectrum function", 

X 1 X -v 

5($, X.) - f f £ f— dvdX . (3) 

1 J 0 X 

If all the Maxwellian orbits were filled (with all velocities and all 
directions) , 5 would be independent of 5 and the profile would be gaussian. 
In this case, however, 5 would also be infinite. In the realistic case, 
wherein dynamical constraints are placed on allowed orbits, the integrated 
density above r^, 

«<»!> * | ~f Kv r , \>«> r . «) 


is finite. Tables exist (Chamberlain 1963, Tables 3 and 4) for computing 

N(X^) for various assumed parameters. 

To evaluate the spectrum function (3) , it is convenient to divide 

particle orbits into three categories: ballistic (gravitationally bound 

and intersecting r ) , satellite (bound but lying wholly above r ) , and 
c c 

escaping. The ballistic plus satellite orbits are collectively called 
"captive", and their integral is rather simple: 


S cap (C ' V 


X -v 

~ (1 - e 2 )dX , 
X^ 


(5) 


where 
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u 2 U> - x- c 


( 6 ) 


Por $ the v integration involve# only the condition of negative 
csp 

energy; it represents the totality of oribts above r, with velocity v f . 

Por ballistic orbits we have the added condition that the orbits 
intersect r c < This integral is more complicated; for reasons that will 
become clear, we carry the functional notation X explicity: 
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where 


and 


v, <0 “ X 2 (X - X + Z 2 ) 
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The totality of all possible satellite oribta is then 
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(9) 


S sat ( V X l } =S cap <*' A l } ~ S bal ( V «' X 1 J ' (l0) 


which yields 


6 


2 + 


S . ( X ; £, X.) 
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- e )dx . (11) 


It must be emphasised that, although X„ does not appear explicitly on 
the right side, when X^ is less than X 2+ (or even X 2 _) , the X^ limit 
prevails. 

The escaping component is analogous to the ballistic; 



for 5 > 0 and 


S e.c a c’ El X l>.’ 0 <13) 

for £ £ 0. Mote that there is a region X 2 _< X < X 2+ that is empty of 

1/2 

escaping particles. However, when £ > X ' /2, both X. and X_, are 

c 2- 2+ 

imaginary and the separate integrals involving these limits in equations 
(7) and (12) merge into a single one and S ** 0. 

Saw 

In practice all possible satellite orbits are no. *r> pected to be 

present. A simple artifice for specifying those present is to define a 

* 

second or "satellite" critical level r , above r . We will assume that 

c c 

any captive particle with a perigee below r * is present; then the 

c 

* 

satellite and ballistic particles are given by S. (X ; £ , X ) , and X 

DAi> OX O 

* 

is changed to X in equations (7) through (11) . 
c 

We now distinguish two cases. First, when the shadow height, r, , 

* 1 

* 

is above r , 
c 


we have simply 
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When the shadow height is belcr t r , we integrate ail bound orbits 

c 

* * 

from r^ up to r c but include above r Q only those orbits intersecting 

* 

r . Thus we write 
c 


SU, *! > l/l - S c>p (S, X x ) - S oat <X c *, E, X c *| 


S «s= ( V 
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(15) 


Equation (14) may be put in the same form as equation (15) by substituting 
%al**c* ; V frotn equation (10) s 


S(E, Xj < X o *, . S cap ({, X x , . S sat ,X c *, E, X x ) 
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* 

Equations (15) and (16) are convenient for computations when X is 

C 

* 

unknown and spectrum functions are to be calculated with varying X to 

C 

obtain a best fit with observations . First S( £ , X^) is computed for 
=0, then S gat may be substracted for any specified . 


III. Sample Profiles 

Throughout the examples shown here, we chose X =7.5, corresponding 

c 

to r being 500 km above the surface of earth and T = 938 for H. The 
c c 

curves in Figure 1 show a range of shadow heights and satellite critical 
* 

levels, r 

c 
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The curves illustrate our contention that good profiles could help 


distinguish the best r as well as the exospheric temperature (s) . For 

c 

* 

the fu) 1 satellite distribution (1 = 0) » there is a drastic cut-off of 

c 

the profile's wings (compared with a gaussian profile) because of the 

absence of orbiting particles exceeding the escape velocity. When the 

* 

satellites are removed altogether (X *7.5), the line center is depleted 
the most, which relatively enhances the wings. The asymmetry of the two 
sides of the profile becomes enhanced, as escaping particles (or their 
absence) become relatively more important. 

Figures 2 and 3 illustrate the two-dimensional phase space that 
contributes to the profiles and may help provide some rationale for why 
the profiles in Figure 1 have such peculiar shapes. The shaded areas 
are forbidded regions — where particles are excluded by the dynamical 
restrictions. These regions are defined by 


£ 2 < (X * - X) 
c 


^ v U ♦ X) 
c c 

. *2 
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(17) 


where 


* , * 2 ., 2 
v = v X /X , 
c c c c 


K < 
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and 
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(18) 


The first restriction requires all orbits to intersect the sphere of radius 
* 

r , whereas the second excludes particles of positive energy from the 
c 

downward-directed hemisphere. 
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legends 

Pig. 1 - Spectrum functions for various shadow heights, r^» and satellite 
critical levels, r , for senith observations from below r. . 

C A 

In all cases the base of the exosphere is specified by X c « 

7.5 (e.g., appropriate for 500 tan altitude on earth, with H 

atoms at 938 \) . For X c * - 1.5 all satellite orbits are 

depleted; for X^ *0, all possible satellite orbits for 

negative total energy are filled, so that particles at any 

1/2 

height are isotropic for all velocities less than X ' . The 
curves, when multiplied by exp (-5?), give the Doppler line 
profiles. 

Fig. 2 - Phase-space diagram showing (by cross-hatched areas) regions in 

* 

which orbits are forbidden for X - X - 7.5. In this caae 

c c 

* 

no satellite orbits are filled; when present (i.e. , when X c 
-> 0) the central forbidden regions are absent, lhe shaded 
areas on the left of each diagram are downward directed orbits 
with velocities exceeding the local escape velocity. When the 
transverse kinetic energy at the critical level, v , exceeds 

v 

the gravitational energy, X , the two forbidden regions overlap. 

C 

On this diagram a particle orbit close to a forbidden region 
is a curve roughly parallel to the boundary of the region. 
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Pig. 3 - Forbidden regions in dimensionless phase space for identical 

transverse kinetic energies (V « 5.0)* A * 7.5* and different 

c c 

satellite critical levels* A * See legend to Figure 2. For 

C 

A^* >7.5 the diagram fits into the sequence of Figure 2. For 

A * < A 2 / 2 v m 5.625, no satellite orbits are forbidden. In 

intermediate cases the forbidden region of satellite orbits 

lies at A<A * and is fixed by the dynamical condition that 
c 

, * 

allowed orbits in this region must intersect A . 
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